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ABSTRACT 

It is well established that Mg II absorption lines detected in background quasar spectra arise from gas struc- 
tures associated with foreground galaxies. The degree to which galaxy evolution is driven by the gas cycling 
through halos is highly uncertain because their gas mass density is poorly constrained. Fitting the Mgll equiv- 
alent width iW) distribution with a Schechter function and applying the A^(Hl)-W correlation of Menard & 
Chelouche, we computed r2(Hl)„j,„ = fi(Hl)haio = l-41!o:44 x 10""^ for 0.4 < z < 1.4. We exclude DLAs from 
our calculations so that il(Hl)haio comprises accreting and/or outflowing halo gas not locked up in cold neu- 
tral clouds. We deduce the cosmic Hi gas mass density fraction in galactic halos traced by Mgll absorp- 
tion is il(Hl)haio/^^(Hl)DLA — 15% and f7(Hl)haio/f^i — 0.3%. Citing several lines of evidence, we propose 
infall/accretion material is sampled by small W whereas outflow/winds are sampled by large W , and find 
17(Hl)infaii is consistent with ri(Hl)outflow for bifurcation at W = 1.23;!;j}28 '-^^'^ accretion would then com- 
prise no more than ^ 7% of of the total Hi mass density. We discuss evidence that (1) the total Hi mass cycling 
through halos remains fairly constant with cosmic time and that the accretion of Hi gas sustains galaxy winds, 
and (2) evolution in the cosmic star formation rate depends primarily on the rate at which cool Hi gas cycles 
through halos. 

Subject headings: galaxies: halos — galaxies: ISM — galaxies: intergalactic medium — ISM: HI — quasars: 
absorption lines 



1. EVITRODUCTION 

Our knowledge of galaxy evolution relies heavily on both 
observations and simulations that focus on the mechanisms by 
which galaxies acquire, chemically enrich, recycle, and expel 
their gaseous component. However, the mean quantity of halo 
gas engaged in any given process remains poorly constrained. 
The use of quasar absorption lines provides a unique tool to 
directly observe these ongoing processes and allows for sensi- 
tive measures of the quantity of gas within galaxy halos along 
with their cross-sections, kinematics, metallicities, densities, 
and temperatures. 

The dense neutral hydrogen that is mostly confined 
within galaxies commonly exhibits damped Lya absorption 
[A^(Hl) > 2 X 10^° cm-2 = DLA]. DLAs are fundamentally 
different from other classes of absorption systems. For exam- 
ple, the Lyman Umit systems [10'^-^^ < A^(Hl) < lO^" -' cm"^ 
= LLS1 have v arying degrees of hydrogen ionization (cf., 
iProchaskall 1 999h and extend out to ^ 100 kpc around galax- 
iefl whereas DLAs are believed to probe the cool, dense pre- 
cursors of star forming molecular clouds and can account for 
up to ^ 50% of the galactic baryonic content (cf., Wolfe et al. 
12001) . 

In an e ffort to de t ermin e the quantity of dense gas within 
galaxies, iRao et all (|200^ measured the cosmological neu- 
tral gas mass density traced by DLAs to be il(Hl)DLA = 
(9.6 ± 4.5) X lO"'* at < z >= 0.92. It has also been 
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argued that fijUl) DLA remains ro ughly constant for z 
0.2-5 (Prochaska & Herbert-Fort 2004; Peroux et al. 2005| 
IRao et al. 2 006; Lah et al. 2008; Son gaila &Cow ie 201 
Meiring et alJboi lb and then decreases by a fa ctor of ~ 2 - 3 
O' ( Z waan et al.l2 005a!; ' Martin et alB " 
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some evide nce showin g that fi(Hl)D 
z = 2.2 (Noterda eme et al.ll2009l) . 

SFR) history I 

has evolved dramatically (e.g.. lBouwens et al.l201 lb and if the 
redshift constancy of r2(Hl)DLA holds, then this implies that 
DLA gas does not directly track the formation of stars. This 
might imply that the global SFR is predominantly governed 
by mechanisms linked to galactic halos, such as gas accretion 
from the intergalactic medium and/or recycling of gas within 
the galaxy halos. An estimate of ri(Hl)haio {excluding Hi 
from DLAs), being the sum of an accreting/infall component 
[r2(Hl)infan] and a wind/outflow component [ri(Hl)outflow], 
could place constraints on the relative importance and roles 
with which these processes drive star formation in galaxies. 

The Mgll AA2796,2803 absorption doublet, which probes 
low-ionization gas with 10'^ < A^(Hl) < 10^^ cm"^, is com- 
monly used to study the gaseous com ponents of galaxies (see 
IChurchill. Kacprzak. & Steidell2005[ for a review). The Mg II 
absorption is observed out to projected galactic radii of ^ 100 
kpc (Kacprzak et al. 2008; Chen et al. 2010a). 

For Mgll absorption systems with rest-frame equivalent 
widths of > 1 A, galaxy color and star for mation rate cor- 
relate s strongly with W (Zibetti et al. 2007; N oterdaeme et alj 
120101: iNestor et alJ |201 L) . a result highly suggestive that 
galactic outflows are responsible for ejecting substantial 
amounts of gas to large galactocentric ra dii. The outflow sce - 
nario is also supported by the result of iBouche et alj (120061) . 
who found an anti-correlation between W and the host halo 
mass, by cross-correlating Mg II absorbers with luminous red 
galaxies from SDSS, and claim that this provides evidence 
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that absorbers are not virialized in gaseous halos of the galax- 
ies. They suggest that the strong absorbers are statistically 
more likely to trace super-winds. 

Indeed, direct evidence for Mg n absorbing winds is seen 
in spectra of star forming galaxies, which exhibit strong out- 
flows blueshifted 300 - 1000 km s~' rela t ive to the galaxy 
jjremonti et al. 200 71: IWeiner et all l2009l iMartin & Bouchd 
l2b09i: .Rubin et al.»2010bl) . These "galaxies almost exclusively 
exhibit W > 1 A absorption. IChelouche & Bowed (120101) 
demonstrated that models of outflowing wind-driven gas re- 
produce the Mg H velocity widths of W > 1 A systems ob- 
served with high resolution. The extent of thes e Mgll absorb- 
ing winds are not well constrained; however, iBordoloi et al.l 
(Hoil) suggest they may reach out to 50 kpc. 

For samp les dominated by W < 1 A, neither IChen et al] 
(l2010a) nor lKacprzak et al.l (1201 Ibl) found a W-g alaxy color 
correlation, contrary to the IZibetti et al.l (120071) result for 
W > I A. Furthermore. Chen et all (1201 Obi) found that Mgn 
"halo size" increases with increasing galaxy stellar mass and 
weakly with specific star formation rate, suggesting a scenario 
in which infalling Mg n absorbing gas structures (selected by 
W < 1 A absorption) fuel s tar formation. 

The SPH simulations of iStewart et al.l (1201 lai) reveal that 
gas-rich mergers and cold-flow streams produce a circum- 
galactic co-rotating, low-ionization gas component that is pre- 
dominately infalling towards the galaxy. In absorption, these 
structures are expected to exhibit ~ 100 km s"' velocity off- 
sets relati ve to the hos t galaxy, consistent wi t h the o bser- 
vations oflS teidel et alj (|2002|), iKac przak et "aL| (|2010a| ), and 
iKacprzak e t al. (2011a). This spatial/kinematic configuration 
yields a correlation between gal axy inclination and W, w hich 
has been previously observed by IKacprzak et al.l (1201 Ibh . 

The above body of evidence suggest that weaker Mgll 
absorption selects gas accretion from infalling cold streams 
or cooled gas returning from earlier processing within the 
galaxy, whereas stronger absorption {W > 1 A) selects out- 
flows. It would be useful to constrain the relative neutral gas 
mass density for both processes, i.e., J7(Hl)haio = ^^(H l)infaii + 
f2(Hl)outflow, to gain insight into how much galactic gas is cy- 
cled through either mechanism at a given time and how this 
compares to J7(Hi)ola- 

In this Letter, we compute the Hi mass den- 
sity within galaxy halos traced by Mgll absorption, 
/^(Hl)haio = rj(Hl)Mi„, . We apply the A^(Hl)-W relation of 
iMenard & Chelouche '(2009) to obtain the Hi column density 
distribution function directly from the Mg II equivalent width 
distribution function. We exclude DLAs from our calcu- 
lations so that i7(Hl)haio budgets gas likely to be accreting 
and/or outflowing from galaxies but not locked up in cold 
neutral clouds. We apply a 1 A bifurcation to compute 
r2(Hl)infaii and ri(Hl)outflow, and determine the W at which 
r2(Hl)infaii = f^(Hl)outflow Throughout we adopt a h = 0.70, 
riivi = 0.3, 51a = 0.7 cosmology. 

2. COMPUTING n(H I) FROM Mg II 

The cosmological neutral hydrogen gas mass density is 
computed from 



17(HI)=^^(A^(HI)).''^ 



C Pc 



dz (l+z)2' 



(1) 



given that. 



(iV(Hl)) — = / Nn(N)dN, (2) 
dz Jo 



where n(N) is the Hi column density distribution function, fi = 
1 .3 is the mean molecular weighQ, Pc is the critical density of 
the universe, Hq is the Hubble constant, m„ is the mass of 
hydrogen, (A^(Hl)) is the mean Hi column density, d'R/dz is 
the number of systems per unit redshift, and E(z) = H{z)/Ho = 
[r!„(l+z)3 + r!A]i/2. 
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Fig. 1. — The di stribution of Mgll res t-frame equivalent w idths, n(W), 
using the data frornlChurchilLetjilJ (i999|) [blue trianglesl. Naraya nan et all 
(2007) [red open squares], Steidel & Sarsent ( 1992) [black open circles] and 
Nestor et al. ( 2005) [green circles]. We fitted the data using a Schechter func- 



tion (solid curve) with parameters C, = 1.08 ±0.12 A 



-0.642 ± 0.062 



and W* = 0.970 it 0.056 A. The dashed curves represent the maximum 1 a 
error in the fit for the parameter uncertainties. 

Computing il{H l) traced by Mg II absorption requires n(N), 
which is not directly observationally known for Mgll ab- 
sorbers. However, Menard & Chelouche (2009) determined 
the geometric mean column density as a function of W to be 
N(lll)=AW^, where A = (3.06 ±0.55) x 10''' cm'^ A"'^ and 
13= 1.73 ±0.26 for 0.5 < W < 3 A and 0.5 < z < 1.4. Em- 
ploying this column density relation, we can obtain n{N)dN 
from the equivalent width distribution, n{W)dW, where n{W) 
is the distribution of Mgll rest-frame equivalent widths. 

In Figure [T] we plot n(W ) versus W us i ng th e 
data o f Steidel & Sargen i (IT992h. iChurchill et alJ (Il999h . 
iNestor et al. (2005), and Naravana n et alJ (l2007f) ~Since the 
work of Steidel & Sargent ( 1992), it has been common to fit 
n{W) as either a power law or an exponential function. For 
W > 0.3 A, .Steidel & Sargent! (119921) c ould not dist i nguish 
which function was preferred, whereas INestor et al.l ( 120051) 
showed that the distribution was clearly exponential. For 

^ 11= 1.3 applies for a fully neutral gas and is the value used for DLA 
studies. For partially ionized gas, fi is slightly smaller, but still of order unity. 
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W < 0.3 A, IChurchill et all (tl999| ) showed that n(W) is a 
power law that is consistent with the lSteidel & Sargeni ( 11992b 
data for < 1 .0 A, which was confirmed by lNarayanan et aP 

§007). 

We used ^-squared minimization to fit a lSchechterl (11976b 
function to the binned data. 



TABLE 1 

f2(Hl) TRACED BY Mgll AT 0.4 < z < L4 



W \ \ W \ / w 
«(W = a( - exp . - 



(3) 



where C» is the normalization such that the unweighted in- 
tegral is equal to d'R/dz (to satisfy Eq. |2]i. In Figure [T] 
we show our fit with parameters C» = 1.08 ±0.12 A~', a = 
-0.642 ± 0.062 and W, = 0.970 ± 0.056 A. The overall equiv- 
alent width distribution at (z) = 1.0 is much better fit with 
a Schechter function than the power law only (see Fig. 6 
of I5lar avanan et al. 2007) or with the exponential only (see 
Fig. 20 of Nestor et al..,2005.) . Interestingly, the characteristic 
equivalent width, W*, marking the transition from a power law 
to an exponential distribution, is consistent with 1 A. Based 
upon the results discussed in ^ we speculate the break in the 
n{W) distribution at W ~ 1 A is physical, that infall/accretion 
structures yield a power law distribution whereas outflow- 
ing/wind structures yield an exponential distribution. 

We note there are subtle differences between the data sets 
shown in Figure [T] The W < 0.3 A data of Churchill et al.j 
(0999) [30 systems] and Narayan an et all (l2007b [112 sys- 
tems] represent 0.4 < z < 1 .4 and have mean redshift of (z) = 
0.9. The W > 0.3 A data of iSteidel & S^genl (11991) [103 
systems] represent 0.2 < z < 2.2 with (z) = 1.1, and the data 
of Nest or et a l. (2005) [1331 systems] representO.4 < z < 2.3, 
also with (z) = 1.1. Thus, the redshift ranges of the low and 
high equivalent width samples are not identi cal, even t hough 
the (z) are fairly consistent. Both Steidel & Sargent! ( 11992b 
and iNestor et al.l (120051) demonstrated redshift evolution of 
n{W) in that there is more power at large W at higher red- 
shifts. This evolution is expected to result in a slightly larger 
W* than if we confined their samples to z < 1.4 (the upper 
limit of the W < 0.3 A data). Given the mean redshifts of both 
of the W > 0.3 A samples, (z) = 1.1, are similar to those of the 
W < 0.3 A samples, (z) = 0.9, we expect that this slightly in- 
creased power in n{W) has only a minor influence on our final 
r2(Hl)haio result, as compared to, for example, the uncertainty 
in the fitted values of A and i n the column density relation 
of iMenard & Chelouch j ( 12009b . We also note (see Figure[T]i 
a slight vertical offset between the two W < 0.3 A data sets 
and between the two W > 0.3 A data sets. In both cases the 
larger surveys yield slightly smaller c/H /dz and the source of 
the offsets is still unknown. 

We now substitute N(Hl)=AW^ into Eqs.[TJ{3]and derive. 



i2(Hl)= -rC^AW^ r(a,w), 

c pc ( l + zY 



(4) 



where r(a,w) is the incomplete gamma function, and where 
a = a + (3+1. The value of w allows integration over specific 
W intervals: w = W,„a.K/W^, for integration w, and w = 
or W,„,„/W* for integration w — > oo, where W„„„ and W„,wc are 
selected cutoffs. 

To compute Eq. |4l we are required to extrap olate the col- 
umn density relation of iMenard & Chelouchd (^009) to W 
values both lower and higher than the domain of their fit. 
Though there is scatter in the column density relation, it cor- 
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" For W > 4.5 A, DLA contamination is 100% and the contribution to f2(Hl)haio 
vanishes. See text for further discussions. 

rectly predicts that W < 0.3 A absorption systems are sub- 
LLS with logMHl) < 17, consistent with the observations of 
IChurchill etan(l2000l) . 

Since we desire to compute the Hi mass density in galaxy 
halos, we correct for the quantity of Hi found in disks. The 
majority of DLAs lie within the optical disk and exhibit cov- 
ering fractions of virtually 100% (e.g., Zwaan et al. 2005b), 
thus correcting for the Hi DLA absorption provides a good 
approximation for the disk contribution of Hi. Though some 
DLAs may arise in halos, this fr action is estimated to be less 
than 1-5% (Fumagall i etani20i 1ft . 

Rao et al. ( 2006) showed that the probability of a Mgll se- 
lected DLA system is P(W) = for W < 0.6 A and then in- 
creases with increasing W for W > 0.6 A. Using a maximum- 
likelihood fit to the binned data in their Fig. 4, we estimate 
this increase as a linear function P(W) = 0.23W- 0.057 for 
0.6 < W < 4.5 A with P(W) = 1 for W > 4.5 A. In or- 
der to correct for DLA contamination in our calculation, we 
weight Eq. H by PiW) - 1. The upper limit of W = 4.5 A 
where DLA contamination is 100%, is consistent with the 
IMenard & Chelouchd (12009b column density relation, which 
predicts logA^(Hl) > 20.3 for W > 4.5 A. 

3. RESULTS AND DISCUSSION 

In Table [T] we present r2(Hl) for selected W ranges. The 
quoted uncertainties are the 1 a confidence levels based upon 
the uncertainties in the fitted parameters C*, a, W*, A, and 
13. We deduce the Hi mass density traced by Mgll absorp- 
tion, which is interpreted as the diffuse Hi contained within 
galaxy halos, is f7(Hl)haio = l-41!o:4l x lO"'^. This value is 
^ 15% of ri(Hl)DLA, indicative that a considerable fraction of 
Hi is contained in galaxy halos relative to the Hi in DLAs; 
it contributes 0.3% to the total baryonic budget (i^i, = 0.045, 
iJarosik et alllMTh at (z) = 1.0. 

We find that W < 0.3 A Mgll absorpti on (often called 
"weak" systems, e.g.. IChurchill et anil999l) selects a small 
fraction of the Hi mass density, r2(Hl)(,<0-3' = 7.7i;|^^ x 10"^ 
From this quantity, it is difficult to ascertain what fraction of 
these systems could be selecting Ly a forest structures versus 
galactic halo structures because estimates of r^Lyc in the ap- 
propriate A^(Hl) range (lO'^'^-lO'^'^ cm"^) are highly uncer- 
tain and are quoted in units of total gas mass density (neutral 4- 
ionize d, cL- lPentoiiet al. 2004). The simulations of lDave elal 
( I20T0I) indicate that n(Hl) ~ 10"^ for A?(Hl) < 10^^ (outside 
halos) at z ^ 0. 1 , where the H I fraction in the Ly a forest is at 
its highest. 

We previously described the theoretical and observational 
evidence supporting the idea that weaker Mg II systems trace 
infall/accretion and stronger systems trace outflow/winds and 
that W* ~ 1 A marks the transition between the two regimes. 
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Applying a 1 A bifurcation to f2(Hl)haio, we find il(Hl)infaii = 
5. 56^1 X 10-5 and r2(Hl)outflow = 8.57!°;^^ x 10"^ The for- 
mer is 6% of fioLA and 0.1% of ilb, and latter is 9% of Hola 
and 0.2% of ilb- The range of W over which the infall and 
outflow J7(Hl) are statistically consistent is W = 1.23;!;Q 2g ^■ 

There is no a priori expectation that our approach to com- 
puting ri(Hl)haio should yield il(Hl)infaii ~ f^(Hl)outflow for 
W ~ W*. Our result may indicate that, over a redshift range 
covering a large percentage of the age of the universe in the 
"post star forming era", a cyclic balance persists between in- 
flow and outflow of galaxies whereby star formation is fueled 
by accreting gas and then an equal mass of gas is ejected back 
into the halos. This is quite suggestive of a halo gas recycling 
m odel (e.g ., Oppenheimer & Dave 2008). 

iNote rdaeme et al. (2009) extrapolated their Hi column den- 
sity distribution below logA^(Hl) = 20.3 and found that LLSs 
contribute c±13% of the total il(Hl) at z > 2.2 (also see 
iPeroux et al.1 llOOSh . If we assume r2(Hl)tot = f2(Hl)haio + 
rj(Hl)DLA, we obtain 17(Hl)haio =^ 13% of rj(Hl)tot for 0.4 < 
z< 1 .4; since f2(Hl)oLA is constant with redshift, this suggests 
that 57(H l)haio has remained constant with redshift and implies 
that the H I mass cycling through halos via infall/outflow has 
also remained constant. 

Thus, given the cosmic evolution of the global SFR (es- 
pecially below z — 2), and presuming galactic infall/outflow 
is strongly coupled to star formation, the global SFR must 
be governed by the rate at which H I gas cycles through ha- 
los (i.e., the SFR and Hi halo cycling rate must evolve in 
parallel). The observation that the mean ionization of Mgll 
absorbers has decreas ed with decreasing redshift from z ^ 2 
( [Bergeron et al.l [T994I) is consistent with this scenario. At 
higher redshift, the more highly ionized H I halo gas selected 
by Mg II absorption constitutes a smaller fraction of the total 
gas associated with Hi; at lower redshift, the more neutral Hi 
gas constitutes a larger fraction of the total gas mass. Thus, 
the total gas mass associated with Hi that is cycling through 
halos is higher at high redshift and lower at low redshift for a 
fixed Hi mass. Based upon simulations, it is predicted that the 
total halo gas mass increases with decreasing redshift to the 
present epoch (e.g., Dave e t al. 19 99), however, this is not in- 
consistent with our proposed scenario because this growth is 
in the "hot" phase that is neither detected via Mg II absorption 
n or a reservoir fo r star formation. 

iRibaudo et al.l (1201 1 ) presented possible observational evi- 
dence of cold accretion in a [Mg/H] = -1.7 LLS at z = 0.27 
near a Z ~ Z© sub-L» galaxy. Cosmological simulations pre- 
dict that cold accretion is truncated at l ow redshifts (e.g., 
iFumagalU et al.l201 UlStewart et al.l201 Ibl) such that the cross 
section of this gas i s a tiny fractio n of the observed Mgll 
cross section (Kacpr zak et al.l 120081; ICheiTe t al. 2010a). If 
cold, metal-poor filaments comprise a component of the in- 
falling material, our findings imply they constitute no more 
than - 7% of l](Hl)tot at (z) = 1 .0. 



Our calculation of ^^( Hl)haio relies heavily on the statisti- 
cal N(Ui)-W relation of iMenard & Chelouchd (120091) . which 
we extrapolated to W = and W = 4.5 A. For our calcula- 
tions of il(Hl)infaii and ri(Hl)outflow, we assumed the break 
in the n{W) Schechter function at W* ~ 1 A is due to in- 
fall/accretion for W < W* and outflowing/winds for W > W, . 
A bifurcation at ~ 1 A in the equivalent width distribution 
separating the two physical processes of infall and outflow 
is an intriguing result that is not well understood within the 
framework of current models. In reality, some fraction of the 
W < 1 A absorbers could arise in winds and some fraction 
of the W > 1 A absorbers could arise in infalling/accreting 
material. Nonetheless, the observational data suggest that the 
majority of W > 1 A are of wind origin, and visa versa. Since 
there are no data to constrain the possible fractional contribu- 
tion of infall or outflow that may reside on either side of W,, 
we make no attempt to quantitatively estimate it. 

Further observations are required to ascertain the veracity 
of this simple scenario. For example, there is mounting evi- 
dence that galax y orientation plays some role in deterniinin g 
W and its origin (iKacprzak et al.l201 IbtlBordoloi et al 1201 lb 
Although, the evidence provided by IKacprzak et al.1 (1201 Ibl) 
supports the idea that W < 1 A systems trace accreting halo 
gas. 



4. CONCLUSION 

We have shown that the Mgll equivalent width distribu- 
tion, n(W), at (z) = 1.0, is well described by a Schechter 
functio n. We combined our n{W) with the .;V(Hl)-W rela- 
tion of iMenard & Chelouchd (12009) to compute ri(Hl) resid- 
ing in galactic halos, as traced by Mgll absorption (excluding 
DLAs). We found that 13% of f2(Hl)tot resides in galaxy ha- 
los and deduced that the infall and outflowing components 
comprise roughly equal Hi mass contributions. The balance 
between the two may suggest that outflows are sustained by 
accretion and that cold accretion by filaments comprises less 
than ^ 7% of r2(Hl)tot. Comparing to high redshift results, it 
appears that f2(Hl)haio has not strongly evolved over cosmic 
time. We argued that this implies that evolution in the cosmic 
SFR must depend primarily on the rate at which cool Hi gas 
cycles through halos, even through the total H I mass cycling 
through halos remains fairly constant. 
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